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Abstract 

While the importance of the seasonal migration of the zonally averaged Hadley circulation 
on interhemispheric transport of trace gases has been recognized, few studies have examined 
the role of the zonally asymmetric monsoonal circulation. This study investigates the role of 
monsoon-like zonally asymmetric heating on interhemispheric transport using a dry atmos- 
pheric model that is forced by idealized Newtonian relaxation to a prescribed radiative equi- 
librium temperature. When only the seasonal cycle of zonally symmetric heating is conside- 
red, the mean age of air in the Southern Hemisphere since last contact with the Northern He- 
misphere midlatitude boundary layer, is much larger than the observations. The introduction 
of monsoon-like zonally asymmetric heating not only reduces the mean age of tropospheric 
air to more realistic values, but also produces an upper-tropospheric cross-equatorial trans- 
port pathway in boreal summer that resembles the transport pathway simulated in the NASA 
Global Modeling Initiative (GMI) Chemistry Transport Model driven with MERRA me- 
teorological fields. These results highlight the monsoon-induced eddy circulation plays an 
important role in the interhemispheric transport of long-lived chemical constituents. 


1 Introduction 


Interhemispheric transport is crucial for the global distribution of long-lived trace ga- 
ses in the troposphere, which are primarily emitted from the Northern Hemisphere (NH) 
midlatitude surface. Cross-equatorial transport has traditionally been quantified by the time 
scale inferred from the two-box model with each hemisphere being a box, using measure- 
ments of rapidly growing, long-lived chemical species such as chloroflourocarbons (CFCs) 
and sulfur hexaflouride (SF¢) [e.g., Heimann and Keeling, 1986; Maiss and Levin, 1994; Le- 
vin and Hesshaimer, 1996; Geller et al., 1997]. The mechanisms of interhemispheric trans- 
port have been interpreted in terms of the seasonal migration of the zonally averaged Had- 
ley circulation [Bowman and Cohen, 1997], the strength and position of tropical convection 
[Lintner et al., 2004], or upper tropospheric westerly ducts associated with Rossby wave bre- 
aking [Staudt et al., 2001]. Recently, Orbe et al. [2016] identified the Asian monsoon anticy- 
clone as an important transport pathway for cross-equatorial transport in the NASA Global 
Modeling Initiative (GMI) Chemistry Transport Model driven with MERRA meteorolo- 
gical fields (referred to as the GMI-MERRA model hereafter). This paper will analyze the 
transport in an idealized atmospheric model in which seasonally-varying zonally symmetric 
heating and the monsoon-like zonally asymmetric heating can be separated, and thus inter- 
hemispheric transport by the seasonally-varying Hadley circulation can be compared with 
that by the monsoonal circulation. A schematic summary of the two mechanisms is depicted 
in Fig. 1, with the boreal winter Hadley cell transporting the NH air mass across the equa- 
tor in the lower level versus the boreal summer monsoon exchanging mass between the two 
hemispheres in the upper level. Since the mass of each hemisphere does not change in the 
time mean, the Hadley cell also transports mass back to the NH in the upper troposphere, and 
the zonally asymmetric monsoon circulation helps maintain the mass balance by either the 
upper tropospheric two-way eddy mixing or the lower tropospheric branch of the monsoon 
circulation. 


Tracer transport in the global atmosphere can be thought of as the consequence of the 
diabatic overturning circulation and isentropic mixing [e.g., Plumb and Mahlman, 1987; Bo- 
wman and Carrie, 2002; Bowman and Erukhimova, 2004; Chen and Plumb, 2014]. In the 
tropics, the probability distribution functions (PDFs) of column water vapor (CWV) and 
trace gases of different chemical sources/sinks exhibit similar exponential tails, providing ob- 
servational evidence of tropical convection’s influence on tracer distributions [Neelin et al., 
2010; Lintner et al., 2011]. Chemical transport models also show that convection can en- 
hance vertical and horizontal tracer transport [e.g., Gilliland and Hartley, 1998; Hess, 2005; 
Erukhimova and Bowman, 2006], and that interhemispheric transport is modulated by the 
seasonal migration of the Hadley circulation [Bowman and Cohen, 1997]. Meanwhile, the 
Asian monsoon anticyclone is associated with horizontal eddy shedding and irreversible 
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(a) Boreal Winter Hadley Cell (b) Boreal Summer Monsoon 
Tropopause Tropopause 
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Figure 1. Schematic summary of the two mechanisms of interhemispheric transport: (a) the boreal winter 
Hadley cell (green) transporting mass across the equator in the lower troposphere, and (b) the boreal summer 
monsoon (red), characterized by the lower-level convergence and upper-level divergence, that exchanges mass 
between the two hemispheres through the upper-tropospheric anticyclone. In both (a) and (b), gray dashed 
lines indicate the annual mean residual meridional circulation, and blue lines depict isentropic surfaces and 


eddy mixing. 


mixing of trace gases along its southward edge [Popovic and Plumb, 2001], and upper tro- 
pospheric Rossby wave breaking is very efficient at mixing air masses across the equator 
[Waugh et al., 1994]. Since tropospheric transport displays transport barriers in the subtro- 
pics [Haynes and Shuckburgh, 2000; Allen and Nakamura, 2001; Chen and Plumb, 2014], 
interhemispheric transport is also quantified by means of a three-box model, in which the 
tropical box is sandwiched between two extratropical boxes [Bowman and Carrie, 2002; Bo- 
wman and Erukhimova, 2004]. 


While the aforementioned advection-diffusion paradigm as well as the two- and three- 
box models are important for our understanding of global tracer transport, the spatial pattern 
of interhemispheric transport is better quantified by the Green function that captures where, 
and after how long, the air in the Southern Hemisphere (SH) was last in the NH. Holzer and 
Boer [2001] examined the mean transport climate, as defined by the Green function with 
respect to localized sources in a climate model, and found that interhemispheric exchange 
times, mixing times, and mean transit times all increase by about 10% under climate war- 
ming. The seasonality of the Green function can also help elucidate the transport pathway by 
which interhemispheric exchange occurs and its relation to both the seasonal movement of 
the Hadley cell and the monsoon [Holzer, 1999]. Furthermore, Holzer [2009a,b] developed a 
rigorous path-density diagnostic that identifies the regions through which air passes from the 
source region (e.g., NH high latitudes) to the receptor region (e.g., SH high latitudes). 


Using surface measurements of SF¢ and CFCs, Waugh et al. [2013] and Holzer and 
Waugh [2015] obtained the mean age and the spectral width of the transit time distribution 
(TTD) for different locations in the troposphere since last contact at the NH midlatitude sur- 
face, which Orbe et al. [2016] then compared with the results from the GMI-MERRA simu- 
lation. They also showed that the TTD in the GMI-MERRA simulation is characterized by 
distinct fast advective transport and slow diffusive recirculations that are deemed important 
for tracers of short and long lifetimes, respectively, and that the upper-tropospheric cross- 
equatorial flow in boreal summer is linked to the Asian summer monsoon. 


The paper is outlined as follows. Section 2 introduces the idealized model of the at- 
mosphere and the diagnostic of tracer transport. The responses in the TTD and its charac- 
teristics to seasonally-varying zonally symmetric heating and the monsoon-like zonally- 
asymmetric heating are compared and discussed in section 3. The transport pathway respon- 
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sible for the simulated TTD is further analyzed in section 4, and the importance of monsoon- 
induced eddy circulations is highlighted. Finally, section 5 provides a brief summary of the 


paper. 


2 Idealized model and transport diagnostics 
2.1 Idealized atmospheric model 


We use the Geophysical Fluid Dynamical Laboratory (GFDL) spectral atmospheric 
dynamical core. The model solves the primitive equations on the sphere with Held and Sua- 
rez [1994] physics: it is forced by Newtonian relaxation to a prescribed radiative equilibrium 
temperature and damped by Rayleigh friction in the planetary boundary layer. Readers are 
referred to Held and Suarez [1994] for details of model physics, and the same Held-Suarez 
parameters are employed here unless otherwise stated. 


We introduce additional perturbations to the Held-Suarez equilibrium temperature to 
mimic the diabatic heating associated with the seasonal cycle in radiative heating and latent 
heat release, as well as the zonally asymmetric monsoon circulation. More specifically, the 
radiative equilibrium temperature is specified as 


oi 
105 
Here ¢ is latitude, p is pressure, x = R/cp = 2/7, where R is the gas constant of dry air 
and c, is specific heat of dry air at constant pressure. The term 20 /sin ¢ gives the zonally 
symmetric off-equatorial heating in summer and cooling in winter. The parameter 7 deter- 
mines the timing throughout the seasonal cycle, which permits to run the model with either 
perpetual equinoctial or seasonally-varying radiative equilibrium temperature. The perpetual 
equinoctial condition is mimicked by setting 7 = 0. To produce a seasonal cycle from winter 
to summer, we set f = cos( ae), where ¢t denotes time. Chen and Plumb [2014] showed that 
the off-equatorial heating can produce a cross-equatorial Hadley circulation and an associa- 
ted wintertime subtropical jet, which can in turn suppress the eddy mixing at the jet core and 
lead to more wave breaking at the jet flanks. 


Teq = max{200, [315 ~ 20 fsin  — 60 sin? @ ~ 10log(=F-) cos? 4]( ar As 


Additionally, 72, specifies a wavenumber 2 perturbation in the NH subtropics to mimic 
the zonally asymmetric effect of Asian and North American monsoons [e.g., Shaw, 2014]. 


Tog = 5 AoH (A) sin(2a) (4. p) “ 


Here H() is a heaviside function. This heaviside function switches on the zonally asymme- 
tric monsoon-like heating in the NH summer but turns it off during the SH summer. The 
perturbation has an amplitude of + Ao with a meridional and vertical structure centered at the 
latitude of ¢9 = 52/36 = 25° and midtroposphere (~600 hPa) 


O(d,p) =. | C6 - dbo)] sine), (Ib dol < , P= 0.2% 10%Pa) 
, 0, (elsewhere) 


Such a meridional pattern has been used to study the impact of zonally symmetric off-equatorial 
heating [Plumb and Hou, 1992] on the cross-equatorial Hadley cell and the influence of zo- 
nally asymmetric heating [Zhai and Boos, 2015] on the monsoon circulation. Here ? = —1 
represents the July condition with the zonally asymmetric heating centered at 25°N. The 
perturbation T;,, is approximately confined in the troposphere by the minimum equilibrium 
temperature of 200K in Eq. (1). The pattern of the radiative equilibrium temperature at the 
summer solstice is displayed for Ag = 60 in Fig. 2. Given the thermal damping time scale 
of 40 days used in the model, this gives a heating rate of ~2 K day~!, which is comparable 
to the magnitude of the zonally asymmetric heating in the observations (Fig. 2 of Hazra and 
Krishnamurthy [2015]). In the boundary layer of the Held-Suarez configuration, the thermal 
damping time scale transitions to very short values near the surface to account for the lack 
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Figure 2. Prescribed radiative equilibrium temperature (K), Teg, at the summer solstice for Ag=60: (a) 
zonally averaged Teg, (b) the horizontal cross section of Teq at 575 hPa, and (c) meridional cross section of 


Teg at 45° longitude. Teq denotes the deviation from the zonal mean. See Eqs. (1) and (2) for details. 
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Table 1. 


Descriptions of idealized experiments. Ag denotes the amplitude of zonally asymmetric heating. 


Pp is the top boundary of the source region Quip. @eddy denotes the strength of eddy mixing. I’ is the mean 


age at 15°S/80°S and 875 hPa since last contact with the NH boundary layer (e.g., 3.1/3.5 denotes the mean 


age of 3.1 years at 15°S and 3.5 years at 80°S). See section 2 for details. 


Expt Brief Description fin Teq Ao in TZ q Pb (hPa) Qeddy IT (year) at 15°S/80°S 
1 no seasonal cycle equinox 0 800 1.0 3.1/3.5 
2a seasonal & zonally symmetric seasonal 0 800 1.0 2.6/3.1 
2b seasonal & monsoon forcing —_ seasonal 10 800 1.0 2.4/2.9 
2c seasonal & monsoon forcing —_ seasonal 20 800 1.0 2.3/2.8 
2d seasonal & monsoon forcing —_ seasonal 30 800 1.0 2.1/2.6 
2e seasonal & monsoon forcing —_ seasonal 40 800 1.0 2.0/2.5 
2f seasonal & monsoon forcing _—_ seasonal 50 800 1.0 1.9/2.4 
2g seasonal & monsoon forcing —_ seasonal 60 800 1.0 1.8/2.3 
3a varied source region depth seasonal 60 700 1.0 1.6/2.1 
3b varied source region depth seasonal 60 900 1.0 2.4/2.9 
4a eddy/mean flow partitioning — seasonal 60 800 0 7.2/6.7 
4b eddy/mean flow partitioning — seasonal 60 800 0.5 2.3/2.5 
4c eddy/mean flow partitioning — seasonal 60 800 2.0 1.6/1.8 


of surface heat and momentum fluxes, which may produce unrealistic diabatic heating rates 
there. 


A series of idealized experiments are conducted and summarized briefly in Table 1. 
Experiment | is simulated under perpetual equinoctial zonally symmetric heating. This is 
in contrast to experiment 2a under seasonally-varying zonally symmetric heating, which is 
used to assess the impact of the seasonal migration of the Hadley cell on interhemispheric 
transport. Next, experiments 2a-2g are performed using the same seasonally-varying zonally 
symmetric heating, but the amplitude of zonally asymmetric heating increases from Ag = 0 
to Ap = 60 at an increment of 10. 


The model simulations are run at T42 horizontal resolution with 20 equally spaced 
sigma levels in the vertical. The subgrid-scale diffusion is parameterized by the V* hyperdif- 
fusion on temperature, vorticity, and divergence, with the diffusion coefficient of 1.0 x 10!° 
m* s~!. For the tracer diagnostics, a grid tracer is used with a semi-Lagrangian scheme for 
horizontal advection [Lin et al., 1994] and a finite-volume parabolic scheme for vertical ad- 
vection. A global mass fixer, that rescales the global tracer mass at every time step [Orbe 
et al., 2012], is used to ensure that the global tracer mass remains unchanged by advection. 
While this does not fix local numerical errors in tracer advection, this ensures global tracer 
mass conservation. There is no explicit numerical diffusion applied to the grid tracer. Note 
that we have used a spectral version of the tracer with explicit hyperdiffusion and find that 
the results are quantitatively similar to the grid tracer. Unlike Orbe et al. [2016], there is no 
convective parameterization in this idealized model and thus no convective tracer transport; 
the short radiative relaxation time scale near the surface may also produce unrealistic near- 
surface diabatic transport. All of the simulations are run for 20 years, with the first 5 years as 
the spinup. 


2.2 Atmospheric circulations 


The atmospheric circulations simulated in the model are exemplified by experiment 
2g in Fig. 3, which is forced by seasonally-varying zonally symmetric heating as well as the 
monsoon-like zonally asymmetric heating (Ag=60). In the NH, the upper-tropospheric zonal 
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Figure 3. Seasonal evolution of atmospheric circulations for the experiment forced by seasonally-varying 


zonally symmetric heating as well as the monsoon-like zonally asymmetric heating for Ag=60 (expt 2g in 


Table 1): (a) anomalous (i.e., annual mean-removed) zonal mean radiative equilibrium temperature Teg (con- 


tours; K) and T/ g at 45° longitude (shading; K) at 575 hPa, (b) zonal mean zonal wind at 375 hPa (m s—!), 


and (c) zonal mean meridional wind at 375 hPa (m s—!), 


mean zonal winds exhibit a seasonal transition from a strong wintertime jet to a weak sum- 
mertime jet, associated with a change from a northward upper-level cross-equatorial flow 

in winter to a southward cross-equatorial flow in summer. The SH displays a similar seaso- 
nal transition with a lag of six months, as expected from the lag in the seasonal cycle of the 
radiative forcing. In addition, the zonally asymmetric heating in the NH subtropics is used 
to mimic the heating anomaly associated with Asian and North American monsoons. This 
yields a hemispheric asymmetry between the summertime circulations in the two hemisphe- 
res, with slightly stronger easterlies and cross-equatorial flow at ~10°N in comparison with 
their SH counterpart at ~10°S. 


The effect of the zonally asymmetric heating is highlighted by contrasting the July 
mean circulations simulated in experiments 2a (Ap=0) and 2g (Ap=60) in Fig. 4. The pres- 
cribed zonally asymmetric heating generates planetary waves that give rise to an upward lift 
in the Hadley circulation near 20°N and enhanced upper-tropospheric angular momentum 
transport from tropical easterlies to midlatitude westerlies. This is, however, secondary as 
compared with the pronounced zonal asymmetry in circulation. The imposed warming leads 
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Figure 4. July mean circulations for the simulations (left) without and (right) with the zonally asymmetric 


heating, i.e., Ag=0 and 60 (expts 2a and 2g), respectively: (top) zonal mean zonal wind (m sy, (middle) 
mean meridional streamfunction (10? kg sy, and (bottom) pressure vertical velocity at 575 hPa (Pa sl), 


Black circles in the bottom right panel indicate the region of local warming (cf. Fig. 2). 


to strong upward motion locally that is accompanied by remote descent to the west of war- 
ming. All these characteristics are consistent with Zhai and Boos [2015], who examine the 
monsoonal circulation response to perpetual local subtropical heating. The streamfunction 
maximum in Fig. 4 is comparable in magnitude with observations (Figs. | and 5 of Dima 
and Wallace [2003]), and the seasonal forcing used here is similar to the sinusoidal temporal 
variation in the leading mode of the Hadley cell. Adding the monsoon-like heating alters the 
structure of the Hadley cell, but it does little to the magnitude of the streamfunction. There- 
fore, it allows us to approximately separate the seasonal migration of the Hadley cell from 
the zonally asymmetric monsoon circulation. However, given the simplicity of the model 
configuration, one should bear in mind that there exist differences in the structure of the Had- 
ley cell between this idealized model and observations. 


2.3 Diagnostic of transport 


While interhemispheric transport has traditionally been expressed in terms of the he- 
mispherically averaged time obtained using a two-box model [e.g., Levin and Hesshaimer, 
1996; Geller et al., 1997], we quantify interhemispheric transport in terms of the mean age 
of air and transit time distribution since last contact at the NH midlatitude surface [Waugh 
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et al., 2013; Holzer and Waugh, 2015; Orbe et al., 2016]. This follows the theory of the age 
of stratospheric air since last contact at the tropical tropopause (see the review by Waugh and 
Hall [2002]), but redefines the boundary condition as the NH midlatitude surface. Unlike 

the path density diagnostic [Holzer, 2009a,b], this does not include the path information of 
where individual air masses pass en route to the SH. 


More specifically, let us consider the transit time distribution corresponding to a con- 
served and passive tracer of mixing ratio y(r, t), where transit time refers to the elapsed time 
Tt =t-’ since the air at a location r and time ¢ was last over the NH midlatitude source re- 
gion at time t’, denoted by Quip and defined here as all gridpoints between 30°N and 50°N 
from the surface to the top boundary at pp = 800 hPa. The tracer mixing ratio in the source 
region Quy at time ¢’ is specified as y(Qmyp, ¢’), and its value is constant in Qyyp to mimic 
the effect of the boundary layer turbulent mixing in more realistic models. There is no source 
or sink of tracer outside the Quip region. Experiments 3a (p,=700 hPa) and 3b (p,=900 
hPa) are conducted to examine the sensitivities of tracer transport to the top boundary of 
Qmaip and are discussed at the end of section 3. 


Given a Green function G(r, t|Qup, ¢’) that satisfies the tracer advection equation and 
the boundary condition G(Qmpp, t/Qmmp, t’) = 6(t — t’), where 6(t — t’) denotes a delta 
function, the tracer mixing ratio at the position r and time ¢ can be written as 


t 
yor?) i (Quan, #”) G(r, tua, 


| X(Qmip, f — T)G(r, T| Quip) dt (4) 


Hence G(r, t| Quo, t’) "propagates" the mixing ratio in the source region Qyyp at time ¢’ to 

the point r at time ft. The Green function or TTD is rewritten as G(r, T}Qurp) = G(r, t]Omnp, f- 
T) under the assumption that the flow is cyclostationary and the dependence on ¢ is drop- 

ped. If the tracer mixing ratio in the source region is specified as y(Qmmp, t’) = 1, it is 
expected that the tracer field will be homogenized eventually, and this yields y(1, 00) = 

ie G(r, T\Quip)dt = 1. Therefore, G(r, t|Qyap)dr can be thought of as the mass fraction 

of the air parcel at the position r and time f that was last in contact with the source region 
QOmnp between the transit time t and t + drt. 


The mean age of air at the position r since last contact with the source region Qyp can 
be defined by the first moment of the TTD 


Tero) =f G(r riQup)dr, (5) 
0 
and the spectral width of the TTD is defined by its second moment 


Acriom) = 4f5 [ [F-TlrlQu)P OC, rlOsnp)ar. 6) 
From a one-dimensional advection-diffusion model with constant advective velocity and dif- 
fusivity, the analytic solution of the TTD can be expressed as an inverse Gaussian distribu- 
tion, G(T) = | = exp[— re} with ¢ = 1t/T, where I, A, and 7 are the mean age, 
spectral width, and transit time, respectively [e.g., Waugh and Hall, 2002]. The TTD of a tra- 
cer subject to pure advection is simplified to a delta function centered at the time scale that 
is determined by the mean advection; enhanced tracer mixing broadens the TTD, leading to 
younger modal age and longer tail of old air (see Fig. 3 of Waugh and Hall [2002]). Hence, 
the ratio A/T, referred to as the shape parameter, can , at least partly, characterize the long 
tail of the TTD in this analytic solution, and thus it will be used to assess the contribution of 
slow eddy-diffusive recirculations relative to fast advective transport. 


The TTD diagnostics are computed in the same way as Orbe et al. [2016] to facilitate 
the comparison with the GMI-MERRA simulation. Idealized Boundary Impulse Response 


(BIR) tracers are implemented to approximate the Green function response to a boundary 
layer delta function forcing. After a 5 year spin-up, four BIR tracers are released in the Quip 
region every three months for a duration of 1 day at the value of | and then held at zero for 
the rest of the integration. Setting the tracer back to zero in the source region acts as a sink 
of the tracer, and thus the BIR tracer decays with time after the initial pulse. For the runs 
subject to the seasonal radiative forcing, four BIR tracers are released on the first day of 
January, April, July, and October, respectively. While the instantaneous variations of TTD 
and BIR are different for a non-stationary flow and the distinctions may be important for the 
annual cycle of the TTD, the TTD and BIR have identical statistics in the long term mean. 
Thus, the average of four BIR tracers centered about the transit time t = 0 is used to estimate 
the annual mean statistics of TTD, although more BIR tracers may be needed to recover the 
exact TTD in a turbulent flow [Haine et al., 2008]. The exponential tail of the TTD near and 
above the tropopause, where the mean age is large, is extrapolated to account for the limited 
length of the simulation. Finally, the accuracy of the BIR approach for computing the mean 
age is verified by a clock tracer, which is held at zero in the source region Quip and is aged 
at arate of 1 year per year outside Quip [Waugh et al., 2013; Orbe et al., 2016]. An ideal 
age is obtained as the 5-year average of the clock tracer after 15 years of integration, when 
the clock tracer has equilibrated almost everywhere in the troposphere. The mean age from 
the BIR tracers agrees quantitatively with the ideal age from the clock tracer, and thus the 
latter is not shown here. 


2.4 Mean mass transport versus eddy mixing 


The effect of mean mass transport versus eddy mixing is diagnosed using the method 
described in Chen and Plumb [2014]. For a tracer of mixing ratio vy, the mean mass trans- 
port and eddy mixing may be separated as 


Ox 


1 1 1 
at = --V- (puy) =--V- (PAmeanUres X) —-V- [PQeddy(U = U-es)X] (7) 
p p p 


where p is density, V is a 3D gradient operator, and u = (u, v, w) is the 3D velocity. For illus- 
tration purposes we assume the flow is incompressible, V - (ou) = 0. Uy-es = (0, Vres, Wres) 
denotes the zonal mean residual meridional circulation, and it is calculated at every time step 
on the model grid to satisfy mass conservation, V - (pU,-es) = 0. As such, the contribution of 
mean mass transport versus eddy mixing to tracer transport can be illuminated by modifying 
the parameters Q@mean and Geddy, respectively. As shown in Chen and Plumb [2014], in the 
limit Of @mean = 1 and @egay = 0, the tracer is advected only by the residual circulation, 
and the tracer is expected to be homogenized along the streamline of the zonal mean residual 
circulation. Conversely, in the limit of @mean = 0 and @aay = 1, the tracer is advected by 
eddy stirring but with no diabatic mass transport, and the tracer is expected to be homogeni- 
zed along isentropic surfaces. By varying the value of @aay from 0 to 2, experiments 4a-4c 
are performed to examine the effect of slow eddy-diffusive recirculations on tracer transport 
in section 4. 


3 Transit Time Distribution (TTD) 


The spatial pattern of interhemispheric transport is described by the TTDs at different 
pressure levels and latitudes, approximated by the average of four BIR tracers released in the 
NH midlatitude boundary layer Qyip on January 1, April 1, July 1, and October 1, respecti- 
vely. Figure 5 gives the TTDs in the upper troposphere (375 hPa) and near the surface (875 
hPa) for experiments 1 (black), 2a (blue), and 2g (red). The latitude 15°S and 15°N are se- 
lected to capture the meridional gradient in the cross-equatorial transport, and 80°S is chosen 
to represent the SH high latitudes. Overall, the TTD in the idealized model is characterized 
by a large peak (modal age) at young transit times and a long flat tail of old air, as found in 
the GMI-MERRA model [Orbe et al., 2016]. In the tropics, the TTDs are heavily skewed to 
young transit times, whereas the TTDs in the SH high latitudes are less skewed and closer to 
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Figure 5. Zonally averaged Transit Time Distributions (TTDs) since last contact with the NH midlatitude 
boundary layer Qyyyp for the simulations (blue) without and (red) with the zonally asymmetric heating, i.e., 
expts 2a (Ag=0) and 2g (Ag=60), respectively, and (black) under the perpetual equinoctial radiative forcing 
(expt 1). The TTD is approximated by the average of four Boundary Impulse Response (BIR) tracers released 


in Quap at the source time ¢’= January 1, April 1, July 1 and October 1, plotted relative to the transit time 


Tt = t — t’, where f is the sample time. The location of each TTD is indicated by pressure and latitude in each 


subplot. Dashed vertical lines indicate the mean age of each TTD. 


the inverse Gaussian distribution with the disappearance of young air (t < 50 days). As such, 
while the modal age captures the peak at young transit times, the mean age is much larger 
than the modal age due to the highly skewed TTD. 


As reviewed in the introduction, interhemispheric transport can be influenced by both 
the seasonal cycle of the Hadley circulation [Bowman and Cohen, 1997] and eddy-diffusive 
recirculations. The former can be attributed to seasonally-varying zonally symmetric hea- 
ting, and the latter may be driven by the monsoonal circulation or some other zonally asym- 
metric forcing. The comparison of different simulations in Fig. 5 shows that the monsoon- 
like zonally asymmetric heating (experiment 2a with Ag = 0 vs. 2g with Ag = 60) is more 
efficient for transport across the equator than the seasonal migration of the Hadley cell (ex- 
periment | vs. 2a). Particularly, the seasonal, zonally symmetric heating reduces the mean 
age at 80°S and 875 hPa from ~3.5 years to ~3.1 years, and the zonally asymmetric heating 
produces additional reduction to ~2.3 years. The larger influence on the age of tropospheric 
air by the zonally asymmetric heating than the seasonal, zonally symmetric heating is also 
evident at other pressure levels and latitudes. 


The spatial structure of TTD characteristics is summarized by four metrics in Fig. 6: 
the mean age I, modal age Tnode, spectral width A, and shape parameter A/T since last con- 
tact with the source region Quyp. In the source region of 30-50°N and below 800 hPa, I = 0 
and A = 0 by definition, that is, the transit time is reset to zero when the isentropic recircu- 
lation travels from the tropics to the NH high latitudes through this region. Hence, in most 
of the NH extratropics this yields a large peak of young air with very short modal age (< 10 
days, noting the logarithmic scale). For the isentropic surface closer to the source region, the 
recirculation is more likely relabeled as young air, and thus the spatial pattern of the mean 
age resembles that of the spectral width, reflecting the change in the long tail of old air with 
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respect to the isentropic surface in the NH. In contrast, in the SH extratropics the gradient 

in the spectral width is small and the spatial pattern of the mean age is more similar to that 
of the modal age, reflecting the gradual disappearance of fast transport paths from the SH 
subtropics to the South Pole (also see Fig. 5). As a result, the spectral width relative to the 
mean age is broad (i.e., A/T > 1) in the NH and narrow (i.e., A/T’ < 1) in the SH. In spite of 
the simplicity of the idealized model, these above features agree well with the GMI-MERRA 
simulation (Fig. 2 of Orbe et al. [2016]) and thus underscore the importance of large-scale 
atmospheric circulations in setting up the transport climate. 


On the other hand, the mean age at the equator displays a sharp meridional gradient 
in contrast to strong vertical gradients in modal age. This vertical gradient in modal age dif- 
fers from the meridional gradient found in the GMI-MERRA simulation [Orbe et al., 2016], 
which can be attributed to the lack of convection and associated vertical mixing in this ide- 
alized dry model that would otherwise mix the modal age vertically by convective transport. 
The similarity of the cross-equatorial meridional gradient between the mean age and spectral 
width in the idealized model again indicates the importance of recirculations in interhemisp- 
heric transport. In more realistic models, interhemispheric transport can be further enhanced 
by convection and vertical mixing [e.g., Gilliland and Hartley, 1998]. 


The TTD response to the monsoon-like zonally asymmetric heating is illuminated 
by comparing the left and right columns of Fig. 6 (experiment 2a with Ap = O vs. 2g with 
Ao = 60). The mean age becomes younger throughout the troposphere, with the largest re- 
duction in the age gradient occurring in the tropics. In contrast to the spatial pattern of cli- 
matological mean age varying with the spectral width only in the NH, this reduction in mean 
age by the zonally asymmetric heating is associated with narrower spectral width throug- 
hout the troposphere. Noting as the amplitude of zonally asymmetric heating increases from 
Ag=0 to Ag=60, the mean age in the SH subtropics/high latitudes at 875 hPa decreases from 
2.6/3.1 years to 1.8/2.3 years, respectively (Table 1). This results in a mean age that agrees 
better with observational-based estimates of mean age of around 1.1 years in the SH subtro- 
pics and 1.4 years in SH high latitudes [Waugh et al., 2013; Holzer and Waugh, 2015] and 
the spatial pattern in the GMI-MERRA simulation [Orbe et al., 2016], although additional 
processes such as convective mixing are needed to produce realistic values in mean age. In 
contrast, there is little change in the modal age and shape parameter except for the modal 
age near 20°N latitude where the zonally asymmetric heating is imposed. These suggest that 
the monsoon-like zonally asymmetric heating has a larger influence on slow recirculations 
(spectral width) compared to fast advective transport (modal age). 


It is worthwhile revisiting the impact of the seasonal migration of the Hadley circu- 
lation on interhemispheric transport [Bowman and Cohen, 1997]. Figure 7 gives the TTD 
characteristics in the simulation in which the radiative equilibrium temperature is fixed to 
the perpetual equinoctial condition, and thus the upward branch of the Hadley cell stays near 
the equator, as evident in the patterns of the mean age, spectral width and shape parameter. 
While the overall patterns of the TTD metrics in the extratropics are qualitatively similar, the 
absence of the seasonal migration in the Hadley cell leads to an older mean age and wider 
spectral width in the SH. This indicates a slowdown in the recirculation of old air, presuma- 
bly resulting from the perpetual updraft at the equator. Again, the seasonal migration in the 
Hadley cell has little impact on the structure of the modal age. 


The robustness of the relationship described above is verified by plotting individual 
metrics as a function of the heating amplitude at selected locations in the tropics and the SH 
high latitudes (Fig. 8). The mean age and spectral width at the selected locations display a 
decreasing trend with a larger heating amplitude, and their decreases in the SH are compa- 
rable in magnitude with the reduction in the NH tropical upper troposphere. In contrast, the 
modal age as well as the shape parameter exhibits little sensitivity to the heating amplitude 
(also see Fig. 5). 
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Figure 6. Zonally averaged mean age (IT; year), modal age (To@e; year), spectral width (A; year), and 


shape parameter (A/T) of the TTDs since last contact with Qyyp for the simulations (left) without and (right) 


with the zonally asymmetric heating (expts 2a with Aj=0 and 2g with Ag=60). Red dotted lines denote the 


280K, 300K and 320K isentropes. Note the log scale in Ty,o¢e- 
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397 Figure 7. As in Fig. 6, but for the simulation under the perpetual equinoctial radiative forcing (expt 1). 
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Figure 8. Mean age (I; year), modal age (Tnode; year), spectral width (A; year), and shape parameter 


(A/T) of the TTDs since last contact with Qygrp as a function of the amplitude of zonally asymmetric heating. 


The location of each metric is indicated in the parentheses as (hPa, degrees latitude). Note the log scale in 
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Figure 9. As in Fig. 5, but for varied top boundaries in Qyrp as pp = 700, 800, and 900 hPa (i.e., expts 3a, 


2g, and 3b), respectively. 


In light of the above role of slow recirculations in determining the mean age and spectral 
width, we further examine the sensitivities of the TTD to the top boundary of the source re- 
gion (Fig. 9). As the isentropic recirculation moves from the tropics to the NH high latitu- 
des, old air is relabeled as young air in the source region, and thus a rise in the top boundary 
of the source region is likely to remove more recirculating old air at the tail of the TTD. In- 
deed, in response to a deeper source region, the TTD exhibits an increase in short transit ti- 
mes and a decrease at the long tail, and the decrease in mean age in the SH is comparable in 
magnitude with that in the NH tropical upper troposphere. Therefore, the similar change in 
the TTD in the SH between Figs. 5 and 9 supports the notion that the increased interhemisp- 
heric exchange in response to the zonally asymmetric heating may be attributed to accelera- 
ted recirculations. 


4 Transport pathways 


The BIR tracers released in four seasons are now compared with one another. This 
allows us to decompose the annual mean BIR or TTD in Fig. 5 by season, and thus to relate 
the TTD change to the seasonal change in the transport pathways associated with the Hadley 
cell or the monsoon [Holzer, 1999; Orbe et al., 2016] 


The seasonal variation of the Hadley cell (Fig. 3) is characterized by the latitudinal 
shift of its upward branch from the SH in January to the NH in July. This latitudinal shift 
influences the location and season of the largest response among the four BIR tracers, as evi- 
dent from the simulation driven by seasonally-varying zonally symmetric heating in Fig. 10: 
the January tracer (i.e., the BIR tracer released at the source time of January 1) exhibits a 
notable increase in the SH tropical lower troposphere (Fig. 10e), whereas the July tracer dis- 
plays a pronounced increase in the NH tropics (Fig. 10c and f), as expected from the latitude 
of the upward branch of the Hadley cell. It is noteworthy that the mean age in the SH tropi- 
cal upper troposphere is smaller than that of the same latitude in the lower troposphere (Figs. 
6 and 7), and thus the primary cross-equatorial transport takes place through the upper tro- 
posphere rather than the lower troposphere. This has been noted by many previous studies 
[e.g., Plumb and Mahlman, 1987; Holzer, 1999; Bowman and Erukhimova, 2004]. It is then 
surprising that the July tracer in Fig. 10e yields the smallest increase among the four seaso- 
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Figure 10. Zonally averaged BIR tracers since last contact with Qyrp for the simulation with seasonally- 


varying zonally symmetric heating (expt 2a). The BIR tracers are released at the source time ¢’ 


(blue), April 1 (cyan), July 1 (red), and October 1 (green), plotted relative to the transit time t = t — t’, where 


t is the sample time. The annual mean is plotted in black. The location of each BIR is indicated by pressure p 


and latitude y in each subplot. 


nal BIR tracers in the SH tropical upper troposphere, in contrast to Orbe et al. [2016] who 
found the July tracer exhibits the largest increase in the SH tropical upper troposphere in the 
GMI-MERRA model (their Fig. 3). 


The inconsistency above can be resolved by the simulation with the monsoon-like zo- 
nally asymmetric heating in Fig. 11. In comparison with Fig. 10e, the January tracer in Fig. 
11e displays a similar, large increase in the SH tropical lower troposphere, but the July tra- 
cer produces a faster response (from ~60 days in Fig. 10c to ~15 days in Fig. 11c) in the 
NH tropical upper troposphere, followed by an enhancement in fast transit times (<60 days) 
in the SH tropical upper troposphere (Fig. 11b). This accelerated upper-tropospheric cross- 
equatorial transport of the July tracer resembles the boreal summer transport pathway found 
in the GMI-MERRA simulation (Fig. 3 of Orbe et al. [2016]). 


The spatial pattern of this July transport pathway is illustrated by plotting the zonally 
averaged July tracer integrated over selected elapsed time intervals since last contact with the 
source region. The left column of Fig. 12 depicts a slow tropical upward transport pathway 
over the time scale of ~1.5 months along the upward branch of the Hadley cell in boreal 
summer. With the aid of the zonally asymmetric heating in the right column of the figure, 
the BIR tracer can reach the upper troposphere at ~20N within 15 days, which, in turn, is 
transported to the SH tropical upper troposphere. 


Figure 13 displays the horizontal distribution of the BIR tracers in the upper and lo- 
wer troposphere integrated over the elapsed time intervals from t =16 to 30 days, when large 
cross-equatorial transport occurs in the zonal mean BIR tracers. In the absence of zonally 
asymmetric heating, the upward transport of the BIR tracer is slow and the subsequent cross- 
equatorial transport is limited. In contrast, the zonally asymmetric heating drives an upward 
motion over the warming center, associated with the convergence in the lower troposphere 
and divergence in the upper troposphere. This up-and-over transport pathway efficiently mo- 
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Figure 11. As in Fig. 10, but for the simulation with the zonally asymmetric heating (expt 2g). 
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Figure 12. Zonally averaged BIR tracers (shading) integrated over selected elapsed time intervals since last 
contact with Qyyp, denoted by [ i BIR(r, t|Qyrp) 4dr, for the simulations (left) without and (right) with the 
zonally asymmetric heating (expts 2a and 2g). The time interval [71, T)] is shown in each subplot. The tracer 
is released at the source time t’ = July 1. Red and blue contours indicate the mean meridional streamfunction 
in July (cf. the middle panel of Fig. 4). 
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Figure 13. BIR tracers (shading) at 375 hPa and 875 hPa integrated over the elapsed time intervals t=[16, 
30] for the simulations (left) without and (right) with the zonally asymmetric heating (expts 2a and 2g). The 


tracer is released at the source time ¢’= July 1, and the horizontal wind vectors represent the climatological 


mean winds in July. Red circles in the right column indicate the region of local warming. 


ves the BIR tracer across the equator through the upper-level anticyclone, which can be assis- 
ted by additional eddy shedding associated with the anticyclone [Popovic and Plumb, 2001]. 


One may ask how much of the mean age of air can be attributed to fast advective trans- 
port versus slow eddy-diffusive recirculations. The eddy mixing in tracer transport can be 
separated from the mean mass transport using a method described in section 2.4. Figure 
14 gives the mean age from the simulations in which the eddy transport is modified by the 
factor @eddy = 0, 0.5, 1.0 and 2.0, and the mean residual meridional circulation is fixed by 
Qmean = 1.0. In spite of the seasonal migration of the Hadley cell, the zero mixing limit le- 
ads to an age distribution as expected from the transport along the mean residual meridional 
circulation. In response to enhanced eddy mixing, the meridional gradient in mean age along 
the NH extratropical isentropic surfaces is reduced; under realistic eddy mixing the isolines 
of mean age become parallel to the isentropic surfaces. The mean age in the tropics remains 
to be vertically aligned irrespective of the strength of eddy mixing. As the largest meridional 
gradient of the mean age is located in the tropics, the change in mean age in the SH follows 
the change in the NH tropical upper troposphere. This confirms the notion that accelerated 
eddy-diffusive recirculations lead to enhanced interhemispheric transport and a reduction in 
mean age in the SH. Interestingly, eddy mixing smooths the separation at ~20°N between the 
vertically-aligned age isoline in the tropics and the isentropic age isoline in the extratropics 
(@eddy=0.5 in experiment 4b vs. @egay=2 in experiment 2g). This is analogous to the impact 
of the zonally asymmetric heating on the mean age distribution at ~20°N in Fig. 6 (Ao = 0 in 
experiment 2a vs. Ag=60 in experiment 2g). 
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Figure 14. Zonally averaged mean age (IT; year) of the TTDs since last contact with Qyrp for the simulati- 
ons in which the strength of eddy mixing is varied by the factor @ggqay = 0,0.5, 1.0 and 2.0 (expts 4a, 4b, 2g, 


4c). See Eq. (7) for details. 


5 Discussion and Conclusions 


While the importance of the seasonal migration of the zonally averaged cross-equatorial 
flow on interhemispheric transport has been recognized [Bowman and Cohen, 1997; Lintner 
et al., 2004], the role of the zonally asymmetric monsoonal circulation has not been well 
understood [Popovic and Plumb, 2001; Orbe et al., 2016]. This study has investigated the 
impact of monsoon-like zonally asymmetric heating on interhemispheric transport, using an 
idealized model of the atmosphere forced by Newtonian relaxation to a prescribed radiative 
equilibrium temperature. The idealized model allows us to separate the zonally symmetric 
heating associated with the seasonal cycle of the Hadley cell from the zonally asymmetric 
heating associated with the monsoon. Interhemispheric transport in the idealized model is 
assessed by the mean age of air since last contact with the NH midlatitude surface [Waugh 
et al., 2013; Orbe et al., 2016]. 


When only zonally symmetric radiative equilibrium temperature is considered, the 
mean age since last contact with the NH midlatitude boundary layer is unrealistically large as 
compared with observational-based estimates [Waugh et al., 2013; Holzer and Waugh, 2015] 
and the GMI-MERRA simulation [Orbe et al., 2016]; the introduction of zonally asymme- 
tric heating reduces the mean age to more realistic values. Interestingly, the modal age dis- 
plays little sensitivity to either the seasonal migration of the Hadley cell or the monsoon- 
like heating, likely due to the lack of convection and vertical mixing in the idealized model. 
Furthermore, the zonal asymmetry produces an upper-tropospheric cross-equatorial trans- 
port pathway in boreal summer that resembles the transport pathway simulated in the GMI- 
MERRA model [Orbe et al., 2016]. This transport pathway is linked to the boreal summer 
monsoon that is characterized by the lower-level convergence and upper-level divergence. A 
distinct pathway associated with the seasonal migration of the Hadley cell, in contrast, trans- 
ports air across the equator to the SH in the lower troposphere in boreal winter. In compari- 
son with Fig. 3 of Orbe et al. [2016], the difference in timing and altitude of the two trans- 
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port pathways (as summarized in Fig. 1) highlights the boreal summer monsoon plays an 
important role in interhemispheric transport. 


The idealized simulations also shed light on the important role of eddy-diffusive reci- 
rculations on interhemispheric transport. The remove of recirculating old air parcels in the 
TTD by raising the top boundary of the NH midlatitude source region leads to a decrease 
in the spectral width and mean age in the SH but with little effect on the modal age. The si- 
milar response in TTD to the monsoon-like zonally asymmetric heating (Fig. 5 vs. Fig. 9) 
suggests that the monsoon may impact the TTD by altering the eddy-diffusive recirculations. 
A direct acceleration of the eddy mixing in tracer transport leads to enhanced interhemis- 
pheric transport and a similar decrease in mean age in the SH (Fig. 14). By analogy, one 
may also argue that convection and associated vertical mixing [Gilliland and Hartley, 1998; 
Hess, 2005; Erukhimova and Bowman, 2006] will ventilate the old air associated with slow 
recirculations and shorten the mean age in the SH since last contact with the NH midlatitude 
boundary layer, a topic warranting future investigation. 
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